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O ■ Abstract 



We present a practical and statistically consistent scheme for actively learning binary 
, classifiers under general loss functions. Our algorithm uses importance weighting to correct 

^ ' sampling bias, and by controlling the variance, we arc able to give rigorous label complexity 

bounds for the learning process. Experiments on passively labeled data show that this 
\l ' approach reduces the label complexity required to achieve good predictive performance on 

many learning problems, 
l/^ \ Keywords: Active learning, importance weighting, sampling bias 



1. Introduction 



' Active learning is typically defined by contrast to the passive model of supervised learning. 

' In passive learning, all the labels for an unlabeled dataset are obtained at once, while in 

active learning the learner interactively chooses which data points to label. The great hope 
\ of active learning is that interaction can substantially reduce the number of labels required, 

\^ • making learning more practical. This hope is known to be valid in certain special cases, 

where the number of labels needed to learn actively has been shown to be logarithmic in 
the usual sample complexity of passive learning; such cases include thresholds on a line, and 
linear separators with a spherically uniform unlabeled data distribution (Dasgupta et al., 
2005). 

Many earlier active learning algorithms, such as (Cohn et al., 1994; Dasgupta et al., 
2005), have problems with data that are not perfectly separable under the given hypothesis 
class. In such cases, they can exhibit a lack of statistical consistency: even with an infinite 
labeling budget, they might not converge to an optimal predictor (see Dasgupta and Hsu 
(2008) for a discussion). 

This problem has recently been addressed in two threads of research. One approach (Bal- 
can et al., 2006; Dasgupta et al., 2008; Hanneke, 2007) constructs learning algorithms that 
explicitly use sample complexity bounds to assess which hypotheses are still "in the run- 
ning" (given the labels seen so far), thereby assessing the relative value of different unlabeled 
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points (in terms of whether they help distinguish between the remaining hypotheses). These 
algorithms have the usual PAC-style convergence guarantees, but they also have rigorous 
label complexity bounds that are in many cases significantly better than the bounds for 
passive supervised learning. However, these algorithms have yet to see practical use. First, 
they are built explicitly for 0-1 loss and are not easily adapted to most other loss functions. 
This is problematic because in many applications, other loss functions are more appropriate 
for describing the problem, or make learning more tractable (as with convex proxy losses 
on linear representations). Second, these algorithms make internal use of generalization 
bounds that are often loose in practice, and they can thus end up requiring far more labels 
than are really necessary. Finally, they typically require an explicit enumeration over the 
hypothesis class (or an e-cover thereof), which is generally computationally intractable. 

The second approach to active learning uses importance weights to correct sampling 
bias (Bach, 2007; Sugiyama, 2006). This approach has only been analyzed in limited set- 
tings. For example, (Bach, 2007) considers linear models and provides an analysis of con- 
sistency in cases where either (i) the model class fits the data perfectly, or (ii) the sampling 
strategy is non-adaptive (that is, the data point queried at time t doesn't depend on the 
sequence of previous queries). The analysis in these works is also asymptotic rather than 
yielding finite label bounds, while minimizing the actual label complexity is of paramount 
importance in active learning. Furthermore, the analysis does not prescribe how to choose 
importance weights, and a poor choice can result in high label complexity. 

Importance-weighted active learning 

We address the problems above with an active learning scheme that provably yields PAC- 
style label complexity guarantees. When presented with an unlabeled point xt, this scheme 
queries its label with a carefully chosen probability pt, taking into account the identity of 
the point and the history of labels seen so far. The points that end up getting labeled are 
then weighted according to the reciprocals of these probabilities (that is, l/pt), in order to 
remove sampling bias. We show (theorem 1) that this simple method guarantees statisti- 
cal consistency: for any distribution and any hypothesis class, active learning eventually 
converges to the optimal hypothesis in the class. 

As in any importance sampling scenario, the biggest challenge is controlling the variance 
of the process. This depends crucially on how the sampling probability pt is chosen. Our 
strategy, roughly, is to make it proportional to the spread of values h{xt), as h ranges over 
the remaining candidate hypotheses (those with good performance on the labeled points 
so far). For this setting of pt, which we call IWAL (loss- weighting), we have two results. 
First, we show (theorem 2) a fallback guarantee that the label complexity is never much 
worse than that of supervised learning. Second, we rigorously analyze the label complexity 
in terms of underlying parameters of the learning problem (theorem 7). Previously, label 
complexity bounds for active learning were only known for 0~1 loss, and were based on the 
disagreement coefficient of the learning problem (Hanneke, 2007). We generalize this notion 
to general loss functions, and analyze label complexity in terms of it. We consider settings 
in which these bounds turn out to be roughly the square root of the sample complexity of 
supervised learning. 
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In addition to these upper bounds, we show a general lower bound on the label com- 
plexity of active learning (theorem 9) that significantly improves the best previous such 
result (Kaariainen, 2006). 

We conduct practical experiments with two IWAL algorithms. The first is a specializa- 
tion of IWAL (loss- weighting) to the case of linear classifiers with convex loss functions; 
here, the algorithm becomes tractable via convex programming (section 7). The sec- 
ond, IWAL(bootstrap), uses a simple bootstrapping scheme that reduces active learning to 
(batch) passive learning without requiring much additional computation (section 7.2). In 
every case, these experiments yield substantial reductions in label complexity compared to 
passive learning, without compromising predictive performance. They suggest that IWAL is 
a practical scheme that can reduce the label complexity of active learning without sacrificing 
the statistical guarantees (like consistency) we take for granted in passive learning. 

Other related work 

The active learning algorithms of Abe and Mamitsuka (1998), based on boosting and bag- 
ging, are similar in spirit to our IWAL (bootstrap) algorithm in section 7.2. But these earlier 
algorithms are not consistent in the presence of adversarial noise: they may never converge 
to the correct solution, even given an infinite label budget. In contrast, IWAL (bootstrap) 
is consistent and satisfies further guarantees (section 2). 

The field of experimental design (Pukelsheim, 2006) emphasizes regression problems in 
which the conditional distribution of the response variable given the predictor variables 
is assumed to lie in a certain class; the goal is to synthesize query points such that the 
resulting least-squares estimator has low variance. In contrast, we are interested in an 
agnostic setting, where no assumptions about the model class being powerful enough to 
represent the ideal solution exist. Moreover, we are not allowed to synthesize queries, but 
merely to choose them from a stream (or pool) of candidate queries provided to us. A 
telling difference between the two models is that in experimental design, it is common to 
query the same point repeatedly, whereas in our setting this would make no sense. 

2. Preliminaries 

Let X be the input space and Y the output space. We consider active learning in the 
streaming setting where at each step t, a learner observes an unlabeled point xt & X and 
has to decide whether to ask for the label yt G Y. The learner works with a hypothesis 
space H = {h : X ^ Z}, where Z is a prediction space. 

The algorithm is evaluated with respect to a given loss function I : Z x Y ^ [0, oo). 
The most common loss function is 0-1 loss, in which Y = Z = { — 1,1} and l{z, y) = l(y ^ 
z) = l{yz < 0). The following examples address the binary case Y = { — 1, 1} with Z gM: 

• Kz,y) = (1 - yz)+ (hinge loss), 

• l{z, y) = ln(l -|- e~y^) (logistic loss), 

• l{z, y) = {y — z)^ = (1 — yz)"^ (squared loss), and 

• l{z,y) = \y — z\ = \1 — yz\ (absolute loss). 
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Notice that all the loss functions mentioned here are of the form l{z,y) = (j){yz) for some 
function (j) on the reals. We specifically highlight this subclass of loss functions when proving 
label complexity bounds. Since these functions are bounded (if Z is), we further assume 
they are normalized to output a value in [0, 1]. 

3. The Importance Weighting Skeleton 

Algorithm 1 describes the basic outline of importance- weighted active learning (IWAL). 
Upon seeing xt, the learner calls a subroutine rejection-threshold (instantiated in later sec- 
tions), which looks at Xt and past history to return the probability pt of requesting yt. 

The algorithm maintains a set of labeled examples seen so far, each with an importance 
weight: if yt ends up being queried, its weight is set to l/pt- 

Algorithm 1 IWAL (subroutine rejection-threshold) 
Set So = 0. 

For t from 1,2,... until the data stream runs out: 

1. Receive xt ■ 

2. Set Pt = rejection-threshold(xf , {xi,yi,pi, Qi : 1 < i < t}). 

3. Flip a coin Qt G {0, 1} with E[Qt] = pt- 

If Qt = 1, request yt and set St = St^i U {{xt, yt, l/pt)}, else St = St^i. 

4. Let ht = argmin/.,g// E(x- .y,c)GSt ^ • l{h{x),y). 

Let D be the underlying probability distribution on X xY. The expected loss of h £ H 
on D is given by L{h) = E(^ ,^)^£) l{h{x),y). Since D is always clear from context, we drop 
it from notation. The importance weighted estimate of the loss at time T is 




where Qt is as defined in the algorithm. It is easy to see that E[Lr(^)] = L{h), with the 
expectation taken over all the random variables involved. Theorem 2 gives large deviation 
bounds for Lt(/i), provided that the probabilities pt are chosen carefully. 

3.1 A safety guarantee for IWAL 

A desirable property for a learning algorithm is consistency: Given an infinite budget 
of unlabeled and labeled examples, does it converge to the best predictor? Some early 
active learning algorithms (Cohn et al., 1994; Dasgupta et al., 2005) do not satisfy this 
baseline guarantee: they have problems if the data cannot be classified perfectly by the 
given hypothesis class. We prove that IWAL algorithms are consistent, as long as pt is 
bounded away from 0. Further, we prove that the label complexity required is within a 
constant factor of supervised learning in the worst case. 
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Theorem 1 For all distributions D, for all finite hypothesis classes H , for any 6 > 0, if 
there is a constant Pmin > such that pt > Pmin for all 1 < t < T, then 



V2 ln\H\+\nl 
max \LT{h) - L{h)\ > ^\ ^ 



< 6. 



Comparing this result to the usual sample complexity bounds in supervised learning (for ex- 
ample, corollary 4.2 of (Langford, 2005)), we see that the label complexity is at most 2/p^;^ 
times that of a supervised algorithm. For simplicity, the bound is given in terms of In \ H\ 
rather than the VC dimension of H. The argument, which is a martingale modification of 
standard results, can be extended to VC spaces. 

Proof Fix the underlying distribution. For a hypothesis h £ H, consider a sequence of 
random variables Ui, . . . , Ut with 



Ut 



Qt 
Pt 



l{h{xt),yt)-L{h). 



Since pt > Pmin, \Ut\ < 1/Pmin- The sequence Zt = X]i=i is a martingale, letting Zq = 0. 
Indeed, for any 1 < t < T, 



E[Zt I Zt^i, 



t,xt,yt,Pt Wt + Zt^i I Zt^i, ... ,Zo] 



Zo\ = 



Pt 



-l{h{xt),yt) - L{h) 



E 



\l{h{xt),yt)-L{h) + Z-, 



t~i I Zt-i 



Zt-i, ■ ■ ■ ,Zq 

■ ■ jZq] = Zt^i 



Observe that \Zt+i - Zt\ = \Ut+i\ < l/pmin for ah < t < T. Using Zt = T{LT{h) - L{h)) 
and applying Azuma's inequality (Azuma, 1967), we see that for any A > 0, 



\LT{h)-L{h)\ > 



A 







= p 


\zr > 




Pmin 



< 2e-^V2. 



Setting A = A/2(ln + ln(2/(5)) and taking a union bound over h £ H then yields the 
desired result. ■ 



4. Setting the Rejection Threshold: Loss Weighting 

Algorithm 2 gives a particular instantiation of the rejection threshold subroutine in IWAL. 
The subroutine maintains an effective hypothesis class Ht^ which is initially all of H and 
then gradually shrinks by setting Ht+i to the subset of Ht whose empirical loss isn't too 
much worse than Lj, the smallest empirical loss in Ht'. 

Ht+i = {h G Ht : Ltih) < L; + At}. 

The allowed slack A* = ^/i8/t) \n{2t{t + 1)\H\^ /6) comes from a standard sample complex- 
ity bound. 
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We will show that, with high probability, any optimal hypothesis h* is always in Ht, 
and thus all other hypotheses can be discarded from consideration. For each xt, the loss- 
weighting scheme looks at the range of predictions on xt made by hypotheses in Ht and sets 
the sampling probability pt to the size of this range. More precisely, 

pt= may: maxl{f{xt),y) - l{g{xt),y). 

f,9&Ht y 

Since the loss values are normalized to lie in [0, 1], we can be sure that pt is also in this 
interval. Next section shows that the resulting IWAL has several desirable properties. 

Algorithm 2 loss-weighting (x, {xi,yi,pi,Qi : i < t}) 

1. Initialize Hq = H. 

2. Update 

= min - — - V— 

heHt^i t - I ^ Pi 



3. Return pt = maxf^g(zHt,yeY Kfi^)^y) - K9{x),y). 



4.1 A generalization bound 

We start with a large deviation bound for each ht output by IWAL (loss- weighting). It is not 
a corollary of theorem 1 because it does not require the sampling probabilities be bounded 
below away from zero. 

Theorem 2 Pick any data distribution D and hypothesis class H , and let h* & H be a 
minimizer of the loss function with respect to D. Pick any 5 > 0. With probability at least 
1 — 5, for any T > 1, 

o h* £ Ht, and 

o L{f) - L{g) < 2Ar_i for any f,g e Ht- 
In particular, if hx is the output of IWAL (loss- weighting), then L[hT) — L{h*) < 21S.t~i- 

We need the following lemma for the proof. 

Lemma 1 For all data distributions D, for all hypothesis classes H , for all 5 > 0, with 
probability at least 1 — 5, for all T and all f,g£ Ht, 

\LT{f)-LT{g)-L{f) + L{g)\<AT. 
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Proof Pick any T and f,g& Hx- Define 

Zt = ^{lif{xt),yt) - l{g{xt),yt)) - {L{f) - L{g)). 

Then E[Zt | Z^, . . . , Zt-i] =^^,,yAlU{xt),yt) - l{9{xt),yt) - {L{f) - L{g)) \ Zi,...,Zt-i] = 
0. Thus Zi, Z2, ... is a martingale difference sequence, and we can use Azuma's inequahty 
to show that its sum is tightly concentrated, if the individual Zt are bounded. 

To check boundedness, observe that since / and g are in Ht, they must also be in 
Hi,H2,...,Ht-i. Thus for ah t < T, pt > \l{f{xt),yt) - l{g{xt),yt)\, whereupon \Zt\ < 
j-\l{f{xt),yt) - Kg{xt),yt)\ + \L{f) - L{g)\ < 2. 

We allow failure probability 6/T{T + 1) at time T. Applying Azuma's inequality, we 
have 



P[|Lt(/) - Lrig) - L{f) + L{g)\ > At] 



f \J2 {^m{Xt),Yt) - l{g{Xt),Yt)) - {L{f) - L{g))^ 



> At 



t=i 



> TAt 



T(r + i)|ff|2' 



Since Ht is a random subset of H, it suffices to take a union bound over all f,g(^H, and 
T. A union bound over T finishes the proof. ■ 



Proof (Theorem 2) Start by assuming that the 1 — 5 probability event of lemma 1 holds. 
We first show by induction that h* = argmin/^g// L{h) is in Ht for all T. It holds at T = 1, 
since Hi = Hq = H. Now suppose it holds at T, and show that it is true at T + 1. Let Ht 
minimize Lt over Ht- By lemma 1, LT{h*) — LT{hT) < L{h*) — L(hT) + A^^ < At- Thus 
LT{h*) < + At and hence h* € Ht+i- 

Since Ht ^ Ht^i, lemma 1 implies that for for any f,g £ Ht, 

L{f) - L{g) < Lt^iU) - LT^i{g) + At-i < L*t_i + At_i - + At-i = 2At_i. 

Since /ir, h* G Ht, we have L(/it) < L{h*) + 2At-i. ■ 



5. Label Complexity 

We showed that the loss of the classifier output by IWAL (loss- weighting) is similar to the 
loss of the classifier chosen passively after seeing all T labels. How many of those T labels 
does the active learner request? 

Dasgupta et al. (2008) studied this question for an active learning scheme under 0-1 
loss. For learning problems with bounded disagreement coefficient (Hanneke, 2007), the 
number of queries was found to be 0{r]T + dlog^ T), where d is the VC dimension of the 
function class, and r/ is the best error rate achievable on the underlying distribution by that 
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function class. We will soon see (section 6) that the term rjT is inevitable for any active 
learning scheme; the remaining term has just a polylogarithmic dependence on T. 

We generalize the disagreement coefficient to arbitrary loss functions and show that, 
under conditions similar to the earlier result, the number of queries is O (^rjT + \/ dT log^ , 
where r] is now the best achievable loss. The inevitable VjT is still there, and the second 
term is still sublinear, though not polylogarithmic as before. 



5.1 Label Complexity: Main Issues 

Suppose the loss function is minimized by h* € H, with L* = L{h*). Theorem 2 shows 
that at time t, the remaining hypotheses Ht include h* and all have losses in the range 
[L*,L* + 2At_i]. We now prove that under suitable conditions, the sampling probability 
Pt has expected value « L* + At_i. Thus the expected total number of labels queried upto 
time T is roughly L*T + Y.J=i ^t-i ~ L*T + ^Tln|i?|. 

To motivate the proof, consider a loss function l{z,y) = (p{yz); all our examples are of 
this form. Say (p is differentiable with < Co < < Ci. Then the sampling probability 
for xt is 

Pt = max max l{f{xt),y) - l{g{xt),y) 

f,g£Htye{-l,+l} 

= max max (t){yf{xt)) - (t){yg{xt)) 
f,9&Ht y 

< Ci max max \yf{xt)-yg{xt)\ 

f,g&Ht y 

= Ci max \f[xt) - g{xt)\ 

< 2Cimax \h(xt) - h* (xt)\. 

heHt 



So Pt is determined by the range of predictions on xt by hypotheses in Hf. Can we bound 
the size of this range, given that any h & Ht has loss at most L* + 2A(_i? 



2At_i > L{h) - L* 

> B^Jlih{x),y)-l{h*{x),y)\-2L 

> B^,yCo\y{h{x)-h*{x))\-2L* 
= CoE^\h{x)-h*{x)\-2L*. 



So we can upperbound maxh^Ht '^x\h{x)—h*{x)\ (in terms of L* and A(„i), whereas we want 
to upperbound the expected value oipt, which is proportional to E^,. max/jg//^ l^ix) — h*{x)\. 
The ratio between these two quantities is related to a fundamental parameter of the learning 
problem, a generalization of the disagreement coefficient (Hanneke, 2007). 

We flesh out this intuition in the remainder of this section. First we describe a broader 
class of loss functions than those considered above (including 0-1 loss, which is not differ- 
entiable); a distance metric on hypotheses, and a generalized disagreement coefficient. We 
then prove that for this broader class, active learning performs better than passive learning 
when the generalized disagreement coefficient is small. 
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5.2 A subclass of loss functions 

We give label complexity upper bounds for a class of loss functions that includes 0-1 loss and 
logistic loss but not hinge loss. Specifically, we require that the loss function has bounded 
slope asymmetry, defined below. 

Recall earlier notation: response space Z, classifier space H = {h : X ^ Z}, and loss 
function I : Z x Y ^ [0, oo). Henceforth, the label space is y = {—1,+!}. 

Definition 3 The slope asymmetry of a loss function I : Z x Y ^ [0, oo) is 

^ ^ max^gy \l{z,y) - l{z',y)\ 

z.z'dZ minygy y) - l{z', y)\ ' 

The slope asymmetry is 1 for 0-1 loss, and oo for hinge loss. For difi'erentiable loss functions 
l{z,y) = (j){yz), it is easily related to bounds on the derivative. 

Lemma 2 Let l^{z,y) = (j){zy), where (p is a differentiahle function defined on Z = 
[—B,B] C M. Suppose Cq < \(j)'{z)\ < Ci for all z £ Z. Then for any z,z' G Z, and 
any y G {-1,+1}, 

Co\z - z'\ < \l^iz,y) - l^{z\y)\ < Ci\z - z'\. 
Thus Iff, has slope asymmetry at most Ci/Cq. 

Proof By the mean value theorem, there is some ^ G Z such that l^(z,y) — l^[z' ,y) = 
(j){yz) - (t){yz') = (t)'{£){yz - yz')- Thus \l^{z,y) - l^{z',y)\ = • \z - z'\, and the rest 

follows from the bounds on (j)' . ■ 

For instance, this immediately applies to logistic loss. 

Corollary 4 Logistic loss l{z,y) = ln(l + e~^^), defined on label space Y = {— 1,+1} and 
response space [—B,B], has slope asymmetry at most 1 + e^. 

5.3 Topologizing the space of classifiers 

We introduce a simple distance function on the space of classifiers. 

Definition 5 For any f,g £ H and distribution D define p{f, g) = E^^£) maxj^ y) — 

K9{x),y)\- For any r > 0, let B{f,r) = {g e H : p{f,g) < r}. 

Suppose L* = mhih^H L{h) is realized at h*. We know that at time t, the remaining 
hypotheses have loss at most L* + 2A(_i. Does this mean they are close to h* in p-distance? 
The ratio between the two can be expressed in terms of the slope asymmetry of the loss. 

Lemma 3 For any distribution D and any loss function with slope asymmetry Ki, we have 
p{h, h*) < Ki{L{h) + L*) for all he H. 
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Proof For any h H, 

p{h,h*) = E^maxy\l{h{x),y) - l{h*{x),y)\ 

< KiE^Jl{h{x),y)-l{h*ix),y)\ 

< Ki {B,,y[l{h{x),y)]+B^,y[l{h*{x),y)]) 
= Ki{L{h) + L{h*)). 



5.4 A generalized disagreement coefficient 

When analyzing the algorithm (Balcan et al., 2006) for active learning under 0-1 loss, 
Hanneke (2007) found that its label complexity could be characterized in terms of what he 
called the disagreement coefficient of the learning problem. We now generalize this notion 
to arbitrary loss functions. 

Definition 6 The disagreement coefficient is the infimum value of 9 such that for all r, 

E3;^DSup;jgB(^. .^)Supy |/(/i(2;),y) - /(/i*(x),y)| < Or. 
Here is a simple example for linear separators. 

Lemma 4 Suppose H consists of linear classifiers {u G M*^ : \\u\\ < B} and the data 
distribution D is uniform over the surface of the unit sphere in W^. Suppose the loss function 
is l{z,y) = (j){yz) for differentiahle (p with Cq < < Ci. Then the disagreement coefficient 
is at most (2Ci/Co)\/d. 

Proof Let h* be the optimal classifier, and h any other classifier with p{h,h*) < r. Let 
u*,u be the corresponding vectors in M'^. Using lemma 2, 

r > E^^DSup \l{h{x),y) - l{h*{x),y)\ 
y 

>CoB^^D\h{x)-h*{x)\ 

= CoEx,~^d\{u — u*) ■ x\ > Co\\u — u* 

Thus for any h G B{h*,r), we have that the corresponding vectors satisfy ||u — < 
2rVd/Co. We can now bound the disagreement coefficient: 

B^^D sup sup\l{h{x),y) - l{h*{x),y)\ 

heB{h*,r) y 

< CiE^^D sup \h{x)-h*{x)\ 

h<^B{h*,r) 

< C7iE^sup{|(n-n*) -xl : \\u-u*\\ < 2r^/Co} 

< Ci ■ 2rVd/Co. 
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5.5 Upper Bound on Label Complexity 

Finally, we give a bound on label complexity for learning problems with bounded disagree- 
ment coefficient and loss functions with bounded slope asymmetry. 

Theorem 7 For all learning problems D and hypothesis spaces H , if the loss function has 
slope asymmetry Ki, and the learning problem has disagreement coefficient 9, then for all 
(5 > 0, with probability at least 1 — 5 over the choice of data, the expected number of labels 
requested by IWAL (loss- weighting) during the first T iterations is at most 

W-Ki- {L*T + 0{^T\n{\H\T/5))), 

where L* is the minimum loss achievable on D by H , and the expectation is over the 
randomness in the selective sampling. 

Proof Suppose h* & H achieves loss L*. Pick any time t. By theorem 2, Ht C {h £ H : 
L{h) < L* + 2Aj_i} and by lemma 3, Ht C B{h*,r) for r = Ki{2L* + 2At_i). Thus, the 
expected value of pt (over the choice of x at time t) is at most 

E^r^D sup sup\l{f{x),y) - l{g{x),y)\ < 2 E^^^, sup sup |/(/i(x), y) - /(/i*(x), y)| 
f,g&Ht y hdHt y 

<2E^^D sup svLY)\l{h{x),y) -l{h*{x),y)\ 
heB{h*,r) y 

< 2er = 46 ■ Ki ■ {L* + At^i) . 
Summing over t = 1, . . . , T, we get the lemma. ■ 



5.6 Other examples of low label complexity 

It is also sometimes possible to achieve substantial label complexity reductions over passive 
learning, even when the slope asymmetry is infinite. 

Example 1 Let the space X be the ball of radius 1 in d dimensions. 

Let the distribution D on X be a point mass at the origin with weight 1-/3 and label 1 
and a point mass at {1,0,0, ... ,0) with weight (3 and label —1 half the time and label for 
the other half the time. 

Let the hypothesis space be linear with weight vectors satisfying \ \w\\ < 1. 

Let the loss of interest be squared loss: l{h{x),y) = {h{x) — y)^ which has infinite slope 
asymmetry. 

Observation 8 For the example above, IWAL(loss-weighting) requires only an expected (3 
fraction of the labeled samples of passive learning to achieve the same loss. 

Proof Passive learning samples from the point mass at the origin a (1 — /5) fraction of 
the time, while active learning only samples from the point mass at (1, 0, 0, . . . , 0) since all 
predictors have the same loss on samples at the origin. 

Since all hypothesis h have the same loss for samples at the origin, only samples not at 
the origin influence the sample complexity. Active learning samples from points not at the 
origin 1/(3 more often than passive learning, implying the theorem. ■ 
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6. A lower bound on label complexity 

(Kiiariainen, 2006) showed that for any hypothesis class H and any rj > e > 0, there is a 
data distribution such that (a) the optimal error rate achievable by H is rj; and (b) any 
active learner that finds h £ H with error rate < r/ + e (with probability > 1/2) must make 
rf' /e^ queries. We now strengthen this lower bound to drf /e^ , where d is the VC dimension 
of H. 

Let's see how this relates to the label complexity rates of the previous section. It is well- 
known that if a supervised learner sees T examples (for any T > d/r]), its final hypothesis 
has error < rj + drilT (Devroye et al., 1996) with high probability. Think of this as r/ + e 
for e = ^Jdrj/T. Our lower bound now implies that an active learner must make at least 
queries. This explains the rjT leading term in all the label complexity bounds 
we have discussed. 

Theorem 9 For any 7],€ > such that 2e < r/ < 1/4, for any input space X and hypothesis 
class H (of functions mapping X into Y = {+1, — l}j of VC dimension 1 < d < oo, there 
is a distribution over X xY such that (a) the best error rate achievable by H is rj; (b) any 
active learner seeking a classifier of error at most r/ + e must make U,{drj^ /e^) queries to 
succeed with probability at least 1/2. 

Proof Pick a set of d points Xo, x\,X2, ■ ■ ■ ^ Xd-i shattered by H. Here is a distribution over 
XxY: point Xq has probability 1— /3, while each of the remaining Xi has probability j3/ (d— 1), 
where (5 = 2(r/ + 2e). At Xq, the response is always y = 1. At Xj, i > 1, the response is y = 1 
with probability 1/2 + 76^, where hi is either +1 or —1, and 7 = 2e//3 = e/(?7 + 2e) < 1/4. 

Nature starts by picking uniformly at random. This defines the target 

hypothesis h*: h*{xo) = 1 and h*{xi) = hi. Its error rate is /3 • (1/2 — 7) = r/. 

Any learner outputs a hypothesis in H and thus implicitly makes guesses at the un- 
derlying hidden bits hi. Unless it correctly determines hi for at least 3/4 of the points 
xi, . . . , the error of its hypothesis will be at least rj + (1/4) • 13 ■ (27) = rj + e. 

Now, suppose the active learner makes < c(d — 1)/7^ queries, where c is a small constant 
(c < 1/125 suffices). We'll show that it fails (outputs a hypothesis with error > r] + e) with 
probability at least 1/2. 

We'll say Xi is heavily queried if the active learner queries it at least 4c/7^ times. At 
most 1/4 of the heavily queried; without loss of generality, these are x\^ . . . , x/^, 

for some k < (d — l)/4. The remaining Xj get so few queries that the learner guesses each 
corresponding bit 5j with probability less than 2/3; this can be derived from Slud's lemma 
(below), which relates the tails of a binomial to that of a normal. 

Let Fi denote the event that the learner gets hi wrong; so Ei^'j > 1/3 for i > k. Since 
k < (ci — l)/4, the probability that the learner fails is given by 

Ppearner fails] = P[Fi + • • • + > {d - l)/4] 

>P[Ffc+i + --- + Fd_i >(d-l)/4] 
> P[B >{d- l)/4] > P[Z > 0] = 1/2, 

where i? is a binomial((3/4)((i — l),l/3) random variable, Z is a standard normal, and 
the last inequality follows from Slud's lemma. Thus the active learner must make at least 
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c{d — l)/7^ = ^{drf /e^) queries to succeed with probability at least 1/2. ■ 



Lemma 5 (Slud (1977)) Let B he a Binomial {n,p) random variable with p < 1/2, and 
let Z be a standard normal. For any k G [np,n(l — p)], P[B > k] > P[Z > {k — 
np)/y/np{l - p)]. 

Theorem 9 uses the same example that is used for lower bounds on supervised sample 
complexity (section 14.4 of (Devroye et al., 1996)), although in that case the lower bound is 
dr}/e^ . The bound for active learning is smaller by a factor of r/ because the active learner 
can avoid making repeated queries to the "heavy" point Xq, whose label is immediately 
obvious. 

7. Implementing IWAL 

IWAL (loss- weighting) can be efficiently implemented in the case where H is the class of 
bounded-length linear separators G M"^ : < B} and the loss function is convex: 

l{z, y) = 4>{yz) for convex cf). 

Each iteration of Algorithm 2 involves solving two optimization problems over a re- 
stricted hypothesis set 

Ht = i^^^h^H■.^Y!il^m^^),y^)<^'+^t']■ 

t'<t 

Replacing each h by its corresponding vector u, this is 

Ht = [^lueM.'^: \\uf < B and ^ ^ ^</)(n • [yiXi]) < L*, + aA . 
t'<t I i=i J 

an intersection of convex constraints. 

The first optimization in Algorithm 2 is = min„g//j, Yli=i ^4'{u ■ {ViXi)), a convex 
program. 

The second optimization is maxtj^„gii/j, (/)(y(u • x)) — (/>(y(v • x)), y£ {+1,-1} (where u,v 
correspond to functions f,g). If (j) is nonincreasing (as it is for 0-1, hinge, or logistic loss), 
then the solution of this problem is max{(j){A(x)) — (f){—A(—x)),(j){A{—x)) — (f){—A{x))}, 
where A{x) is the solution of a convex program: A{x) = miiiu^HT ^ " ^- The two cases 
inside the max correspond to the choices y = 1 and y = —1. 

Thus Algorithm 2 can be efficiently implemented for nonincreasing convex loss functions 
and bounded-length linear separators. In our experiments, we use a simpler implementation. 
For the first problem (determining L^), we minimize over H rather than Hj-', for the second 
(determining A(x)), instead of defining Ht by T — 1 convex constraints, we simply enforce 
the last of these constraints (corresponding to time T — 1). This may lead to an overly 
conservative choice of pt, but by theorem 1, the consistency of /ly is assured. 
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7.1 Experiments 

Recent consistent active learning algorithms (Balcan et al., 2006; Dasgupta et al., 2008) have 
suffered from computational intractability. This section shows that importance weighted 
active learning is practical. 

We implemented IWAL with loss- weighting for linear separators under logistic loss. As 
outlined above, the algorithm involves two convex optimizations as subroutines. These were 
coded using log-barrier methods (section 11.2 of (Boyd and Vandenberghe, 2004)). We tried 
out the algorithm on the MNIST data set of handwritten digits by picking out the 3's and 
5's as two classes, and choosing 1000 exemplars of each for training and another 1000 of 
each for testing. We used PCA to reduce the dimension from 784 to 25. The algorithm 
uses a generalization bound of the form ^ d/t; since this is believed to often be loose 
in high dimensions, we also tried a more optimistic bound of In either case, active 

learning achieved very similar performance (in terms of test error or test logistic loss) to a 
supervised learner that saw all the labels. The active learner asked for less than 1/3 of the 
labels. 
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Figure 1: Top: Test logistic loss as number of points seen grows from to 2000 (solid: 
supervised; dotted: active learning). Bottom: ^queries vs ^points seen. 



7.2 Bootstrap instantiation of IWAL 

This section reports another practical implementation of IWAL, using a simple bootstrap- 
ping scheme to compute the rejection threshold. A set H of predictors is trained on some ini- 
tial set of labeled examples and serves as an approximation of the version space. Given a new 
unlabeled example x, the sampling probability is set topmin+(l— Pmin) \'^^'^;hi,h.j&H L{hi{x),y) — 
L{hj{x),y)] , where pmin is a lower bound on the sampling probability. 

We implemented this scheme for binary and multiclass classification loss, using 10 de- 
cision trees bootstrapped on the initial 1/lOth of the training set, setting Pmin = 0.1. For 
simplicity, we did't retrain the predictors for each new queried point, i.e., the predictors were 
trained once on the initial sample. The final predictor is trained on the collected importance- 
weighted training set, and tested on the test set. The Costing technique (Zadrozny et al., 
2003) was used to remove the importance weights using rejection sampling. (The same tech- 
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nique can be applied to any loss function.) The resulting unweighted classification problem 
was then solved using a decision tree learner (J48). On the same MNIST dataset as in 
section 7.1, the scheme performed just as well as passive learning, using only 65.6% of the 
labels (see Figure 2). 
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Figure 2: Top: Test error as number of points seen grows from 200 (the size of the ini- 
tial batch, where active learning queries every label) to 2000 (solid: supervised; 
dotted: active learning). Bottom: ^^queries vs ^points seen. 



The following table reports additional experiments performed on standard benchmark 
datasets, bootstrapped on the initial 10%. 



Data set 


IWAL 


Passive 


Queried 


Train/test 




error rate 


error rate 




split 


adult 


14.1% 


14.5% 


40% 


4000/2000 


letter 


13.8% 


13.0% 


75.0% 


14000/6000 


pima 


23.3% 


26.4% 


67.6% 


538/230 


spambase 


9.0% 


8.9% 


44.2% 


3221/1380 


yeast 


28.8% 


28.6% 


82.2% 


1000/500 



8. Conclusion 

The IWAL algorithms and analysis presented here remove many reasonable objections to the 
deployment of active learning. IWAL satisfies the same convergence guarantee as common 
supervised learning algorithms, it can take advantage of standard algorithms (section 7.2), 
it can deal with very flexible losses, and in theory and practice it can yield substantial label 
complexity improvements. 

Empirically, in every experiment we have tried, IWAL has substantially reduced the 
label complexity compared to supervised learning, with no sacrifice in performance on the 
same number of unlabeled examples. Since IWAL explicitly accounts for sample selection 
bias, we can be sure that these experiments are valid for use in constructing new datasets. 
This implies another subtle advantage: because the sampling bias is known, it is possi- 
ble to hypothesize and check the performance of IWAL algorithms on datasets drawn by 
IWAL. This potential for self-tuning off-policy evaluation is extremely useful when labels 
are expensive. 
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